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In this review, we will discuss how the chiral symmetry and Ua(1) breaking effects are 
reflected in the correlation functions. Using the Banks-Casher formula, one can identify 
the density of zero eigenvalues to be the common ingredient that governs the chiral 
symmetry breaking in correlation functions between currents composed of light quarks 
with or without a heavy quark. Similarly the presence of the 1/4(1) breaking effect is 
determined through the contribution of the topologically non trivial configurations that 
depends on the number of flavors. We also discuss how the symmetry breaking effects 
are reflected in the gluon correlation functions. Finally, we review the Witten Veneziano 
formula for the n' mass in medium. 



1. Introduction 

The breaking and its possible restoration of a symmetry at finite temperature 
and/or density is a fascinating subject as the possible restoration in medium could 
be probed in relativistic heavy ion collision and/or in nuclear target experiments. 
The restoration of chiral symmetry has been linked to the vector meson spectral 
density an d has been the subject of great theoretical and experimental interest up 
to this dajffl^El n has also been linked to quenching of the pion decay constant!^ 
and possible observation of the sigma meson in nuclear matter through the n — tt 
correlation^^. As for the Ua(V) symmetry, its breaking is due to the anomaly, which 
induces an operator relation that remains broken above the QCD phase transition. 
However, the physical effect, such as the large rf mass, is intricately related to chiral 
symmetry breaking and the question of whether the mass will remain constant near 
the chiral symmetry restoration point is of particul ar intere st as the partial quench- 
ing could be observed in nuclear target experiment [P I ^ HUI Moreover, the two pion 
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Bose-Einstein correlation o bserved at RHIC seems to suggest the quenching of the 
r)' mass at high temperatur JHH2H31 

In this work, we will deliver a unifying picture on how the symmetries are 
reflected on the correlation functions. In particular, we show through the Banks- 
Casher formula" that one can identify the density of zero eigenvalues to be the 
common factor that governs t he c hiral symmetry in all order parameters constructed 
from the correlation functions^! Similarly the presence of the Ua{X) breaking effect 
is determined through the contribution o f the t opologically non trivial configura- 
tions that depends on the number of flavorPEZ] We also discu ss h ow the symmetry 
breaking effects are reflected in the gluon correlation functions^! Furthermore, we 
show how these results can be used to generalize the Witten Veneziano formula to 
predict the properties of r]' mass in mediunPSl. 

2. Chiral symmetry and zero eigenvalues 
2.1. Banks- Casher formula 

The expectation value of an operator in QCD is defined as follows: 

(0(x)) = i I W]0{x) (1) 



z . 

where the positive definite measure is defined as 

Z = J[dA„{x)]exp[—^ J d 4 xF 2 }det[l?+m} = J [da]. (2) 

The notations and definition of the Euclidean path integral are introduced in the 
appendix. 

It has been shown that from a Euclidean path integral point of view, the origin of 
chiral symmetry breaking is the p rese nce of zero eigenvalues for the Dirac equation 
in the presence of the gauge fields^! 

(qq(0)) = -<Tr[(0|— 3— 10)]> = -7r(Tr[J A=0 (0, 0)]> 
L/> + rn 

= -tt< J ^Tt[J x=0 (x, x)}) = -n(p(X = 0)), (3) 
where we define a current density matrix of zero eigenvalues as follows: 

J\=o{x,y) = Y,Mx)1>l(.v)- ( 4 ) 

Here, ilftipx = <^A> an d the sum is over zero eigenvalues only. p(X) is the density of 
eigenvalue A. Relation Eq. (J3J shows that the essence of chiral symmetry breaking 
is the non- vanishing zero mode current densit y matr ix Ja=o(2:, y), which reduces to 
the scalar density of zero eig envalues p(X = 0ji4H5| m Eq. ©. 
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2.2. Correlation functions of light meson 

2.2.1. Connected diagrams 

There are other operators and/or correlation functions that reflect the symmetry 
structure of the vaccuum/medium. One can show that the presence of non- vanishing 
density of zero eigenvalues is the unifying ingredient that dictates the breaking 
of chiral symmetry in these operators. To show this, let us study the differences 
between meson two point functions of chiral partners. The simplest example is the 
difference between scalar-pscudoscalar two point functions. 



, s .,(q) = / d A xe i ^(rUr a q(x)qr a q(0) - qr^qix^i^q^J) 

d A xe lqx (Tr [r Q r Q (i-f 5 S(x, 0)i"f 5 S(0, x) - S(x, 0)5(0, x))]) 



d 4 xe iqx (Ti 



-T 2 (0| 



1 



\x) 



Lp + m lp + m 



d 4 xe iqx (Ti 



r 2 (0|- 



1 



\x) x J o (x,0) 2ir 



(5) 



Therefore, the difference vanishes when the current density matrix of zero modes 
Jo vanishes as in Eq. ((3|). It should be noted that the argument of the wave function 
are identical in Eq. while it is not so in Eq. (JS|), the difference in x will introduce 
a correlation length that is not of interest at this point. 

The difference between the vector and axial vector also has the common factor 
(S(x,y) — i'j 5 S(x,y)ij 5 ) given in Eq. (|A.21j) . which is proportional to the density 
matrix of zero modes. Eq. (|A.21|) itself is gauge variant, but becomes gauge invariant 
in the correlation function. 



^v-A 



(?) 



d 4 xe lqx {T[ qT a ^q(x)qr a ^q(0) - qr a 7 5 ^q(x)qr a 7 5 ^q(0) 



d 4 xe lqx (Tr[r 2 j 5 ^S(x, 0) 7 5 7m S(0, x) - ^S(x, 0) 7m S(0, x))}) 



m— >0 



d 4 xe lqx (Tr 



T 7„(0|tt— k)7^o(*,0) 2tt 
lb + m 



(6) 



As with Eq. ([5]), Eq. ([5]) is non- vanishing when Jo ^ and so are all order param- 
eters of chiral symmetry breaking. 
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2.2.2. Disconnected diagrams 

If there is no flavor matrix in the scalar part of Eq. ([5]), such as in the SU(2) case, 
disconnected diagrams will contribute to the correlation function. 

As-p(9) = J^xe^(T^qq(x)qq(0)-^i^q{x)qr a i^q(0)^) 
d 4 xe iqx {Tr[S(x,x)]Tr[S(0,0)]) + connected terms 



d 4 xe lqx (Tr 
/^<Tr 



ttJo(x, x) 



\x) 



Tr 



Tr 



(°lnT^I°) 



ty + m 

7T J (0,0) 



+ 



Hence, the disconnected part also depends on J$^. 



(7) 



2.3. Correlation functions of Baryons 

For the case of a nucleon, a typical interpolating field couples to both the positive 
and negative parity nucleon. S everal methods have been developed and attempted 
to isolate the parity eigenstate o 1QI2QI. The problem becomes simple, if one is inter- 
ested in the differences in parity partners of the nucleons. Here, one can look at 
few examples of differences between currents of chiral partners of the nucleon. For 
example, an interpolating current for A could be n^ — e abc (u^ij 5 Cdi,)h c ; where h 
is a heavy quark. Then, we have 



*A-A 



(?) 



d 4 xe^e abc e a ' h ' c V^^7 5 C'4)^(a;)(S a ^ 7 5 C ,T 4)/i c ,(0) 

(ulCd c )h c (x)(u a ,CdJ,)h c ,(0) 

d 4 xe lqx 6(Tr[S d (x, 0)C T ^(x, 0)C - S d (x, 0)z 75 C T ^(x, 0)Cij 5 ]S h (x, 0)) 
1 



d 4 xe iqx 6(Ti 



(x\ 



d 4 xe iqx l2ir(Tr 



-\0)[(x 
S^0,x)Mx,0) 



— |0)-z 7 5 (x|— ^|0)i 7 5 



S h (x,0)), 



S h (x,0)) 
(8) 



where the chiral limit is taken with the light quark mass. Again, one notes that the 
difference is proportional to Jo. 

When one replaces the heavy quark h to the u quark in the first two lines in 
Eq. ([8]) as in the case of the proton correlation function, there is an additional 
contraction contributing to the correlation function. This extra term is given as 
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follows: 
A™ P *(<Z) 



d 4 xe* qx 6(S u (x,0)((x 



m->-Q 



IJ)+m 

d 4 xe lqx 12ir(S u {x, 0)J Q (x, 0)S u (x, 0)), 



|0)-*7 5 W- 



-I/>+m 



|0)i 7 a )S u (x,0)) 



(9) 



which is again proportional to Jo. 

One can choose another form for the current tja = ( abc (u^Ch^i^ 5 d c . Then, the 
difference between chiral partners takes the following form: 



A 



A— A* 



(?) 



^ xe ^^ e «^^(rK<C , ft 6 )*7 6 dc(!c)(«a'^)dc'(0)*7 8 

(ulCh b )d c {x){u a ,Chl)d c ,{Q) 

d 4 xe iqx 6(Tr[S h (x, i))C T Sl{x, 0)C] ( S d (x, 0) - i~f 5 C T Sj(x, 0)Cij 5 



S h (x,0)S*(0,x) 



J o (x,0)), 



(10) 



which is again proportional to Jo- In this case, the current with h replaced by u 
would be identically zero as the current itself vanishes. 



2.4. Correlation functions of meson with one heavy quark 

Chiral symmetry is also reflected well when there are heavy quarks. Consider the 
meson correlation function in Eq. ([5]). If one of the quarks is heavy, the correlation 
function becomes 

A 5-p0?) = / d 4 xe lqx (T[hq(x)qh(0) - hi-f 5 q(x)qij 5 h{0))) 



(11) 



d 4 xe lqx (Tr[ ( ij 5 S(x, 0)if 5 S h (0, x) - S(x, 0)S h (0, x))]) 
1 



d 4 xe lqx (Ti 



m h 



\x)J o (x,0) 



Again, the difference vanishes when the current density matrix of zero modes Jo 
vanishes as in Eq. 

If one further takes the heavy quark limit nih — >• oo of Eq. (|11[) . One finds 



A s _p(?) m =°^°° I d 4 xe lqx (Tr 
= — <W>. 



(0| 



Ip + m h 



\X)J (X,0)[27T 



} 



(12) 



One can obtain Weinberg type sum rules^ for the heavy light quark system 
by saturating the spectral representation of Eq. (fTT|) with the ground states of the 
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corresponding channels and then equating it with Eq. (I12[) . Then, assuming that 
the contributions from the excited states and/or continuum cancel, one finds 

fD m D ~ fD m D = 2m h{qq), (13) 

where the /d,/d and mD,mD are the decay constants and masses of the pseu- 
doscalar and scalar heavy light meson. The expected scaling of the decay constants 
and masses with the heavy quark mass are^l 

t -1/2 

f D oc m h , 

mo cx mh- (14) 

As expected, Eq. (fT3j) is consistent with the heavy quark scaling, as the light quark 
condensate on the right hand side does not scale with the heavy quark mass. More 
sum rules involving higher moments can be obtained by expanding the last line of 
Eq. (|11[) with respect to the heavy quark mass. 

We can use Eq. (|13p for the charm sector to estimate fr> . Using fn = 206.7 
MeV, m D = 1865 MeV, m Do = 2318 MeV^, TOc = 1.26 GeV and (qq) = 
(0.225GeV) 3 used in QCD sum rulei^, one finds that 



f Do = 182 MeV. 



(15) 



This is slightly smaller than fjj, as was found to be so in a QCD sum rule 
cal culatio n^] Lattice calculation also seems to point to a slightly lower value for 



1261271 



. Similar sum rule can be obtained for the vector and axial vector me- 



son and their corresponding decay constants. The generalization to medium of the 
formula in Eq. (|13p will te nd to modify the masses and the decay constants when 

JOQI 

chiral symmetry is reduced^ 21 . 



2.5. Correlation functions of Baryon with one heavy quark 

In the heavy quark limit — > oo, Eq. ([8]) simplifies to the following: 



A A-A* (?) 



m fe-y o 127T 



St (0,0) Jo (0,0) 



Noting that 

ReS* (0,0) = ReS(0,0) = vrJ o (0,0), 
the real part of Eq. (IT6|) simplifies to the following: 

.oo 12tt 2 



CI(«) 



■(Tr 



Jo(0,0)J o (0,0) 



(16) 
(17) 
(18) 



It should be noted that the operator form of Eq. (|T5)) can be obtained directly from 
the first line of Eq. (JSJ . Equating this one finds 



127r 2 (Tr 



Jo (0,0) Jo (0,0) 



-{u T a Cd c ){u al C^)) 



(19) 



I 



I 
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2.6. I^-quark operators 

One also notes from a direct calculation that left hand side of Eq. (fT9|) is equal to 
a four quark operator composed of quark- ant iquark pairs as given in the following 
identity. 

= 2((u a ij 5 d a )(d b i^f 5 u b ) - {u a d a ){d b u b )) (20) 

Hence, one can obtain identity constraints among chiral operators. 

The identity of the four quark operators appearing in Eq. (|20[) follows from the 
fact that both are the same chiral order parameters. There are other four quark 
operators that are related to the chiral order parameters. For example, the four 
quark operator appearing in the vector and axial vector meson sum rule can be 



combined into a four quark operator that is also a chiral order parameter^!. 



((qT 1 X a T b q)(qr 2 \ a T b q) - {qT^ X a T b q)(qi 1 b T 2 X a r h q)) 



-(Tr 



AVr2 ( '^' ) r ^ '^'°)-^( '^' ^ 



AVr 2 S(0, 0^^(0,0) 2vr 



> 

)• (21) 



Here, A a can either be the color matrix and/or the unit matrix. r a is the flavor 
matrix. The inclusion of the flavor matrix is important to exclude the disconnected 
diagram. It should be noted that if the color and flavor matrix are normalized as 
A 2 = 1, then the matrix element is independent of whether we include the color 
matrix or take it to be just the unit matrix. Also, Fi, T 2 can be any gamma matrix 
with either identical and/or different Lorentz index. It should be noted that Eq. (|21l) 
can also be written in terms of the left and right handed quarks as follows: 

2((qr 1 X a r b q R )(q L r 2 X a r b q) + (qr 1 X a T b q L )(q R T 2 X a T b q)). (22) 

Several identities follow from Eq. ([21]) or Eq. ([22]) . For example, 

((uX a d)(dX a u) - (uij 5 X a d)(di-f 5 X a u)) = ((ud)(du) - (uij 5 d)(dij 5 u)) 

= 2((u L d R ){d L u R ) + (d R u L )(u R d L )), 

(23) 

which vanishes when chiral symmetry is restored. One can take Ti and r2 to hav e 
Lorentz indices and derive identities involving twist-4 operator J^^^^ I ^2 | 33 | 34 | 35 | 

3. Ua(X) effect and Topological configurations 

The integration over gauge fields given in Eq.Q can have nontrivial configurations 
characterized by v = (g 2 /327r 2 ) J d 4 xFF ^ 0. The QCD partition function can be 
divided and summed over by different topological configurations. 

Z = Y J Z V . (24) 
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Configuration with nontrivial topological charges have n+ {n—) number of right 
handed (left handed) zero mode solutions that satisfy the index theorem v — 
n + — ri-. For these configurations it is useful to explicitly write the zero mode 
contribution with asymmetric chirality as a separate quark determinant. For exam- 
ple, for configurations with v = 1, the partition function is given as follows: 

Z„=i = j [dA^(x)] v =i[dq(x)dq(x)]'[dq L dq R ] exp - S Y m - j d A xq{x)[I]> + m]q{x) 

= j [dA ll (x)] 1/= i exp[-SyM]dct'[i/> + m] x dct( J d i x^p (x)rmpo(x)^ 

= j '[d/i] 1/= idet( J d i x-ip (x)rmf; a (x)^ . (25) 



Here, tpo is the zero mode solution and the measure in the zero mode runs over all 
flavors. The prime after the quark measure means that the zero mode is taken out; 
the integration over these quarks give the primed quark determinant. The second 
determinant comes from the integral over the zero modes. One should note that 
the contribution from these configurations is zero in the chiral limit. Therefore, the 
path integral for the difference in the correlation functions discussed in the previous 
section can be thought of the integration over gauge field configurations with v = 0. 

However when calculating the n-point functions involving quark operators with 
asymmetric chirality containing all flavors, the topological configuration will give 
non vanishing contributions. This is so because in such £/a(1) variant configurations, 
the integral over the chiral zero modes are saturated by the external operators and 
do not appear in the zero mode part of the quark determinant. 

Let us consider the 2 flavor (u,d quarks) case. For the quark condensate, 



(qq(0)) = (qq(0)) v = + [dA fl (x)] l/=1 [dq(x)dq(x)}' 
x [du L du R ][dd L dd R ] I q L qR(0) I exp 



Sym - J d 4 xq(x){IJ> + m]q(x) 
Y j [4*W2(^oM0))(y" d 4 xMx)mMx?) + (26) 



As can be seen in the second term of the above equation, the topological configura- 
tions are proportional to the light quark mass and do not contribute to the quark 
condensate. However, when one looks at the two point function, the zero mode 
integral can be saturated by the external operators: 



(qq(x)qq(0)) = {qq(x)qq(0)) l , =o 
1 
Z 

z 



^\ f [ d v]v=i[ du R}[ du L}i dd R\[ dd L] (qLqii{x)j ^<7L<?i?(o)^ + • • 

+ 4 / [d/iUi2 (i>oMx)) (Mo(0)) +■■■ (27) 
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The topological configurations will therefore contribute when looking at the differ- 
ence between correlation functions related by the Ua(\) transformation. For exam- 
ple, 

(qq(x)qq(0)) - (qT 3 q(x)q T 3 q(0)) 

= (~)v=a + \ J [M^ykMx)j UoMO)) + ■■■ (28) 

For the case of Nf flavors, v = 1 configuration will still be proportional to 
m Nf ~ 2 and not contribute to the four quark operators. In this case, the lowest 
dimensional operators where the v = 1 configuration contribute are the 27V f quark 

. JTH1 

operators^. 

Summarizing the sections so far, one can conclude that chiral symmetry breaking 
occurs when the current density matrix of zero modes are non zero, while the Ua(1) 
breaking effect comes about when the topological configuration contributes to the 
operator expectation value and thus in principle independent of chiral symmetry 
breaking. Phcnomenologically, it means that even if chiral symmetry is restored 
at finite temperature or density, the /7a(1) breaking effect will persist in rt-point 
functions and will lead to phcnomenologically observable consequences. However, as 
discussed before, how big "n" has to be depends on the number of flavor, and how 
exact the statements are depends of the magnitude of explicit symmetry breaking. 
One of the obvious physical observable is the mass of a hadron. Therefore, we will 
look at how the relations on the J7a(1) effect is related to the rf mass. 

4. Gluonic correlation 

4.1. LET for Gluonic operators 

Let us first review the derivation of the low energy theorems for gluonic 



correlational. We define the scalar and pseudo-scalar gluonic two-point function, 
S(q) = i J Ae" 1 (t^ G 2 (x)^G 2 (0)^ , 

P(q) = i [ d A xe l(lx ( T—GG(x)—GG(0)\ . (29) 

J \ 4lT 47T / 

Here, G 2 = G^G^ and GG = G^^e^apG^p- There is a well known low energy 
theorem that follows from noting that the measure given in Eq. ([2]) leads to the 
following identity: 

.(O) =i [ d 4 x(TO(x)g 2 G 2 (Q)), (30) 



d(-i/V) 

where the bare coupling g 2 is related to the ultraviolet cutoff Mo to the unique 
scale of QCD. 

A = M exp(~V (31) 



i 
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where 6=11 — |iV/ 133. If there are no other scales in the system, the physically 
observable matrix element should be proportional to the d'th power of the scale, 
where d is the dimension of the operator. Therefore, 



d(-W) 



^^rfn7v)( constxAd ) 



32tt 2 

Substituting O — (3a s /47r)G 2 to the left hand side of the first line of Eq. (|2l)]) leads 
to the following low energy theorem for the scalar correlator. 

TT^^TT^^G 2 ) = * I d 4 x(T^G 2 (x)^G 2 (0)}, (33) 

167T 2 G?(-l/4# 2 ) 47T J 47T An 

The situation for the pseudo-scalar is more involved. The low energy theorem 
depends on whether the correlation function is calculated in pure gauge theory or in 
the presence of light quarks. When there are light quarks present, the pseudo-scalar 
gluon current can be written as 



N f a s 
Air 



Therefore, for Nj — 3, P(q) can be written in terms of the quark currents. 



P(q) =q»q v i J d 4 x e^ x (Tq^ l5 q(x) g 7 ,75<z(0)). (35) 

Hence, 

P(q = 0) = (36) 

However, when there are no light quarks, there are similar low energy theorems as 
in the scalar gluonic correlator. To derive a similar equation for the pseudo-scalar 
current, first introduce a 9 term as in Eq. (|C.1[) . Then we can have expectation 
value of parity odd operators. In such a theory, 

(^GG)e = i J d 4 x(r^GG(x).g 2 G 2 (0)) 9 , (37) 



tf(-l/4 5 2 ) x 4tt 

Using the heavy quark mass expansion, we can write it into the following form, 

<%{-%)m h hv%h) = i f d 4 x(r^GG(x)(4 7 r 2 )(-12)m, i M(0)) ( ,(38) 



d(-l/4. 9 2 )V ' "'° J v 4tt 

As explained in Appendix C, we make a chiral rotation of the heavy quark by 7r/2: 
This will only change the heavy quark mass matrix to rrih — > 771^175. Then Eq. (|38l) 
will change to the following: 

{-mm h hh) 6 = 1 I d 4 x(T^GG(x)(4n 2 )(-12)m h hi l5 h(0))o.(39) 



d(-l/4 5 2 )M v ; n '° J Air 



1 



1 
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Making use of the heavy quark mass expansion again, we obtain the following 
relation for the low energy theorem for the pseudo scalar current. 

1 d ,3a, . f , A , 3a s „ 3a, 



12tt 2 rf(-l/4 9 2 ) V 4tt 



■G 2 ) = i J d*x(T^GG(x)^GG(0)), (40) 



where we have taken the 9 to be zero at the end. Although we have used the heavy 
quark mass expansion, this result should be valid as we can take the mass to be 
infinitely heavy. Using the formula in Eq. (|32[) . we finally obtain the low energy 
theorem for both the scalar and pseudo scalar gluonic correlators. 

S( q = 0) = ^G 2 ) 

P (<Z = 0) = -f(^G 2 ), (41) 

7T 

where the subscript in the pseudo-scalar current represents calculation in pure 
gauge theory. 

At finite temperature, the constant gets modified because there are extra scales 
in the system; temperature and/or density. Then the matrix element could depend 
directly on the extra scale. Hence, one has to make sure that the derivative in 
Eq. (|30|) do not act on these part. Such operations can be taken into account 
explicitly by noting that the matrix element can now be expressed as follows: 

(0) U . T = const X A d f(j,^j, (42) 

where the subscript a, T denote that the expectation value is taken at finite chemical 
potential and/or temperature. Eq. (|30|) will then be modified as 

and the low energy theorem will look as follows: 

p ^=v=-l( d - T w->%K G2 > m 



4.2. Gluonic correlator in medium 

Now let us go back to the pseudo scalar correlation function given in Eq. (|2T)1) and 
discuss their fate when chiral symmetry is restored. Here, we will again introduce 
light quarks with Nf flavors. Then using Eq. (|34p. the pseudoscalar gluonic current 
can be written in terms of the divergence of the axial current. 



P(q) =q»q v i \ d 4 xe lq x 



(Tqjaj5q(x) 97^75^(0)) - (Tq^q(x) qi v q{0)) 



(45) 



i 
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where we have subtracted out the zero contribution from the conserved vector 
current. Using the previous representations, the correlation function can be written 
as 



P(q) =q»q" J d 4 xe lq - x 



(Tr[S(x, x) 7M 7 5 ]Tr[S(0, 0) 7 „7 5 ] - ^[S(x, x)^]Ti[S(0, 0) 7 ,]> 



- (Tr[S(0, x)^ l5 S(x, 0)7,75] + Tr[5(0, x)^S(x, 0) 7 ,]} 



(46) 



Now, when chiral symmetry is restored, the two terms in the second line of the 
above equation will cancel each other, as they are the same as the difference between 
flavored chiral partners; vector and axial vector currents. The remaining first line 
constitutes the disconnected contributions. However, as discussed before and in 
appendix Eq. (IA.22|) . the disconnected pieces all vanish in the chiral limit when 
chiral symmetry is restored. This is again because the disconnected contributions 
are proportional to the current density matrix of zero eigenvalues. 

Tr[S A (x,x)} ~ Tr[5 A (x,x)r] = Tr[J A=0 (x, x)T] -> 0, (47) 

where T is a Hermitian gamma matrix-!^. Therefore, when chiral symmetry is re- 
stored, 

P(q) -> 0, (48) 

in the chiral limit for any finite external momenta q. It should be noted that when 
Eq. (|48p is true only when chiral symmetry is restored, while Eq. (I36[) is always 
true. 



5. r\' mass in medium 

We have discussed the relationship between chiral symmetry and XJ a(X) effects in 
both the vacuum and in medium. When chiral symmetry is restored, the correlation 
functions between chiral partners will become degenerate. In addition, the two 
point correlation functions that are related by V and chiral transformations 
will become degenerate when the number of massless flavors is larger than two. For 
the physical case, this means that the effect of ?/a(1) symmetry in the two point 
functions will only be proportional to the strange quark mass. Here, we will review 
how such relation can be related to the fate of the T)' mass when chiral symmetry 
is restored^. 

5.1. Witten- Veneziano formula 

One can dire ctly re late the pseudo-scalar current and the mass using the large 
N c arguments ^ 39 !. 
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In terms of the physical states, the pseudo-scalar gluonic correlation function 
looks as follows: 

|(Q||gGG|n th glueball)| 2 
P® = - E jTTMl 



E 



| (0 1 If GG | n th meson)| 2 
q 2 - ml 



= P (q) + Px{q). (49) 

The first term in the right hand side of the above equation indicates contributions 
from glueballs, while the second term shows those from the mesons composed of light 
quarks. As discussed before, P(q = 0) = while Po(q = 0) ^ when we assume 
massless quark. This seems in contradiction to the large N c argument because the 
Pq(q = 0) that scales as N 2 is canceled by quark effects that scales as order N c . It 
was noted that this cancelation is possible by the existence of the rj' whose mass 
scales as l/N c and cancels the gluonic effect in Eq. ([49]) only when q = 0; other 
meson masses have a smooth large N c limit. This constraint directly relates the r/' 
mass to the low energy theorem. 

P 0(0 ) = JWfWI'. (50) 



By using the J/a(1) anomaly, 

,„, 3a 



y GG\r,') = (0\d M Jg\rf) (5i) 



Eq. ([50]) becomes as follows: 

Po(0) = -mfrfiNf, (52) 

where Nf is the number of light flavors. In Eq. ([ST]) , we made use of /,,' = f n to 
lowest order in 7V C . Eq. (f5"2"]) is the celebrated WV formula. 
Substituting Eq. (glj) into Eq. ([52]), we find 



"V = \ ^U-G 2 ) 1/2 « 464 MeV, (53) 
V 33 U I" 

where we have used f n = 130 MeV, (^G 2 ) = (0.35GeV) 4 and used b = 11 for pure 
glue theory. This is smaller than the vacuum value of the rj' mass as expected; this 
is the part coming from the 11^(1) effect to the mass of 77'. 



5.2. Witten- Veneziano formula in medium 

To obtain the generalized formula at finite temperature, one notes that the thermal 
gluonic effects are of order iV 2 while that of the quarks are of order N c . If one is in 
the confined phase, the phase is composed of mesons, glueballs and baryons: The 



1 



1 
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scattering of these states scale as order 1, 1 and N c respectively. Hence as long as one 
assumes that the number of hadrons do not scale with N c , hadronic effects can be 
neglected. Near the phase transition, the degeneracy of hadrons would increase and 
follow the scaling of gluons and quarks. Therefore, the leading order effect would 
come from the gluons. Taking into account the effect of the thermal gluons, the 
low energy theorem given in Eq. (|44[) will be modified only by the effect of finite 
temperature. Therefore, the generalization of the WV formula given in Eq. (|50|) 
would be the following. 

ml, = t[ — ~, 4^ ■ 54 

" rid- T-^) (^G 2 )„ V ' 

b \ al J \ 7T I T,puro gauge 

In obtaining the result, we have assumed that the general structure given in Eq. (|49l) 
does not change. Naively, finite temperature effect will introduce the medium rest 
frame and break Lorentz invariance; this means that the correlation function is 
a separate function of qo and q. However, the thermal correction we consider is 
embedded in the scalar gluon condensate and in the scalar low energy theorem. 
Therefore, the corresponding change should also be only in the scalar quantity; 
hence Eq. (|49|) should be valid for our purpose. 

First let us consider the denominator of Eq. (|54[) . It has been known for a long 
time, that the gluon condensate has contribution from the perturbative and non- 
perturbative parts. Moreover, it was also known that at the critical temperature, 
the non-perturbative contribution changes abruptly, but does not y anish completely, 
and retains more than half of its non-perturbative value I 40 | 41 | 42 ] 

The effect of subtracting out the second term in the denominator of Eq. (|54| 
is to get rid of the perturbative correction, or the seemingly scale breaking effect 
that is not related to scale breaking but due to the introduction of an external 
scale parameter T. T he lea ding perturbative correction to the gluon condensate is 
proportional to gf T 4 1 43 | 44 | x nere f orej assuming that the temperature dependence 
is of the following form, 



we find, 



—G ) =G (T) + agtT\ (55) 

7T 



' l - T wf){ : f G % = [ d - T w^ a ' {T] 



if the temperature dependence of g s is neglected. The only temperature dependence 
that survives is Go(T), whose scale dependence is coming from dimensional trans- 
mutation and not from the external temperature only. It is the non-perturbative 
part that dominates the behavior of the right hand side of Eq. (|54|) . Moreover, as we 
discussed before, for the gluon condensate, one has to use the lattice result obtained 
in the pure-gauge theory calculation. The critical temperature in such a calculation 
Tp U re-gauge ~ 260 MeV is known to be around 100 MeV larger than that from a 
full QCD calculation Tqcd^- On the other hand, while the expected change of the 



i 



i 
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gluon condensate is more abrupt in the pure-gauge calculation, the actual change 
in the condensate value itself at the critical temperature is found to be similar to 
the full calculation^!. This means that the change of the gluon condensate can be 
effectively neglected up to temperatures near Tqcd- 

Finally, we use the fact that when chiral symmetry is restored, Eq. (|54p is zero 
for large q as given by Eq. (|48|) . Hence, when chiral symmetry is restored, 

{Q\GG\rf) ~ 0{m q ). (57) 

Therefore, going back to Eq. ([54)) and making use of the previous discussions, we 
find that when chiral symmetry is restored, 

m n , — > 0, (58) 

in the chiral limit. One concludes that in the large N c limit of QCD, rj' mass 
will become degenerate with the other Goldstone bosons. A similar conclusion was 
obtained in Ref. SSI that the anomalous Ua(1) rf mass squared vanishes at high T 
as the chiral quark condensate (qq) . 



6. Conclusions 

We have used Euclidean path integral to show that the current density matrix of 
zero modes is responsible for the breaking of chiral symmetry in operators con- 
structed from correlation functions. This is a generalization of the Banks-Casher 
formula to any order parameter. As for the ?7a(1) symmetry, the contribution of 
topologically non trivial configurations determine the presence of its breaking ef- 
fect. The contribution of topological configurations however depend on the number 
of flavors and the number of external legs of the correlator. We have also discussed 
how the mechanism can be applied to the gluon correlation functions. Combining 
the results, we have described how the Wittcn Vcneziano formula for the rf mass 
can be generalized to the medium. 
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Appendix A. Notations 

Here we will summarize the notations. The QCD partition function is given by 

Z = J [dA fl (x)}[dq(x)dq(x)] exp[iS QC D}, (A.l) 



where 



Sqcd = J d A x 



jF 2 + qW(d tl -iA ll )-m q )q 



(A.2) 



i 
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The Wick rotation to imaginary time, 

xq = t — > —ir = —ixi (A. 3) 

comes with the following: 

d A x -> -id A x E (A.4) 

A -> iA A (A. 5) 

V" = -Fl + F% -> Fl + 4- (A.6) 

-> -7o(9 4 - iAi) - hi(di - iAi) = - iA^lA.7) 

where = 70 , — i-fi . It should be noted that the gauge field here is redefined 
such that the field strength tensor is independent of the coupling. One can go 
back to the usual definition of the gauge field and field strength tensor G M „ by 
redefining the following substitution. 

The partition function then becomes 



Z = J [dA fi (x)][dq(x)dq(x)} exp[-S% C D], (A.9) 

where 

\w + <H;„- 

and = — igA^. Also, Sym will be used for the gluonic part only. 

After integrating over the quark field, the partition function in Eq. ()A.9|) can 
be written as, 



S QCD = I d4x E 



(A.10) 



Z = J [dA^(x)} exp[-S yM ]dct[# + m] 

[df,] (A.ll) 



(1) One should note that d/i is a positive definite measure. To prove this, one notes 
that the eigenvalues of the determinant comes in pairs. If ip is an eigenvector 
iljhjj = Xtp, then so is 7 5 V> because il/>j 5 ip = —ij^Ifop = —Xj 5 tp. Hence, 

det[# + m] = Y[{m - iX) = JJ (m 2 + A 2 ) > (A.12) 

A A>0 

(2) Sometimes it is useful to write the determinant part as 

Z = J[dA ll (x)]dek[l? + m\ cxp[ — Sym] 

[dA^x)]e T ^^ m » exp[-S YM ] (A.13) 



i 



i 
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(0(x)) = \ I [da]0(x) (A.14) 



(3) The vacuum expectation of an operator 0(x) is defined as 

1 

~Z 

A.l. Derivation of Banks- Casher formula 

Here, we repeat the derivation of the Banks-Casher formula ^31 that is extensively 
used in numerical simulations to calculate chiral symmetry breaking. The chiral 
order parameter is defined as follows: 

(qq) =yj d?x(qq(x)) =~J d^Tr^)^))] 
= -i/^(Tr[(x|^|,)]> 

= -\\ [d M ]Tr[(0|— j— 10)], (A.15) 
Z J IJ> + m 

where in the first line, the measure given in Eq. (IA.9[) has been used. One notices 
that the quark propagator selects to divide out the particular matrix element from 
the quark determinant. 

To further calculate the trace of the inverse quark propagator, one can use the 
solutions of the Dirac equation and express the chiral condensate in terms of the 
wave functions as follows: 

ijj\(x)^ ijj\(x) 1 



m = - f d ^ Mx 7 x{x) 1 - J 

J t—- 1 V m — iX 

= -<E^ A )^-T>- (A.16) 

To discuss chiral symmetry breaking in the chiral limit, we have to look at 
lim m _j.o (<7<z) • To this end, we note, limm^o m -i\ = ^(t) + ttS(X). Now, the 
imaginary part is zero because the eigenvalues are paired (positive and nega- 
tive eigenvalues) and the eigenfunction of the pairs are related by 7 s so that 
^_ x (x)i()-\(x) = ip\(x)ip\(x). Hence, 

(qq) = -n(J2^{mx(0)) = -7r(Tr[J A=0 (0, 0)])>. (A.17) 

A=0 

where the current density matrix is defined in Eq. (j4]). A simple booking way to 
understand the meaning of X)a=o * s ^° assume that the original sum of eigenvalues 
consisted of continuous eigenvalues and a degenerate set of finite zero eigenvalues. 
The remaining sum is the sum over finite set of remaining zero eigenvalues. This 
relation shows that the essence of chiral symmetry breaking is the non-vanishing 
current density matrix of zero eigenvalues. For a system of finite box of volume V, 



1 



1 
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the zero modes are small eigenvalues that vanish in the infinite box limit. In practice, 
the wave functions of zero modes are expected to be localized around a topological 
gauge configuration and hence chiral symmetry breaking can be translated to non 
vanishing expectation value of the matrix ?/>o(O)'0o(O) ■ 

Another derivation comes from using Eq. (IA.9l) and Eq. (IA.10[) so that, 



V 



( / d 4 xqq(x)) 



1 dZ 



V Zdm a 



1 

V 



d 4 xTr[(a 



(A.18) 



where the last line follows using Eq. (|A.13J) . 

The disappearance of the imaginary part in Eq. (|A.16j) can be implemented at 
the operator level by the following method: 



(qq) = -(Tr 



(0| 



1 



Ip+m 
1 



10) 



> 



(0| 



Ip + m 



|0)-* 7 5 (0| 



^VA(0)Vi(0)(2vr ( 5(A) 



- > 

2 y 



= -7r(Tr[J A = (0,0)]). 

The identity that we will be using extensively is 
1 



Ip + m 



Ip + m 



\y) 



(A.19) 



(A.20) 



and 



(x\- 



-\y) - h 5 (x\- 



m— >0 



^VA(x)Vi(2/)(2vr ( 5(A) 



J\=a{x,y) \ 2?r 



(A.21) 



This relation itself is not gauge invariant, but still provides the essential interme- 
diate step that relates chiral symmetry breaking to correlation functions. 



A. 2. Chiral symmetry breaking 

As can be seen in Eq. (jA~2T|) . the essence of chiral symmetry breaking is the non 
vanishing current density matrix of zero eigenvalues J\=o(x, y). That is when chiral 
symmetry is restored, any expectation value proportional to the density matrix will 



i 



i 
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vanish. For example, the disconnected diagram of a two point function is propor- 
tional to the following trace: 

MzlraT— Wr] = ^[J x=0 (x,x)r], (A.22) 

Lp -j- 777 

where T is any gamma matrix coming from the current of the two point function 



A. 2.1. Gluon condensate and eigenvalues 

Some insights can be obtained by looking at the heavy quark condensate, which in 
the heavy quark mass expansion can be related to the gluon condensate. 

(hh) = — (tt^-G 2 ) (A.23) 

777/1 1Z7T 

(hh) = ~ ( d 4 x(Tr[(x|— L-lar)]) 
V J IJ> + m h 

= - [ d *x(Y Mx) ' Mx) 1 ) 

J X V V m-iX' 

-)• (A.24) 



777 7 ^ 



A 

ip\(x)^ip\(x), 



Therefore, 



;^G 2 Hl2<]Tp(A)}. (A.25) 

7T » — ' 



Appendix B. Topological configurations 

The integration over gauge fields given in Eq.([5]) can have nontrivial gauge config- 
urations characterized by v — (g 2 /32ir 2 ) J d 4 xFF ^ 0. These configurations have 
77 + (77-) number of right handed (left handed) zero mode solutions that satisfy 
the index theorem v = n + — 72_. For these configurations it is useful to explicitly 
write the zero mode contribution as a separate quark determinant. For example, 
for configurations with v = 1, 

Z v= i = J [dA^xJlv^i exp[-SYM][dq(x)dq(x)][dq R dq L ] exp[- J d A xq{x)[If> + m]q(x)} 



= j [dA,j,(x)] v=1 cxp[- / —^F 2 )det'[tf> + m} x det( / d i xip (x)mip (x) 



1 

2? 

= J [d/i]„=idet( J d i x-ijjo(x)mipQ(x)\ (B.l) 

Here, tpo is the zero mode solution and the measure in the zero mode runs over 
all flavors, ipo is the zero mode solution. One should note that the contribution 
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from these configurations is zero in the chiral limit. However calculating the n- 
point functions involving quark operators with asymmetric chirality containing all 
flavors, the quark determinant from the zero modes does not appear and topological 
configuration will give non vanishing contributions. 



Appendix C. Adding the 9 term 

Consider adding the 6 term to the QCD action. 

Z g = [dA fl (x)][dq(x)dq{x)] 



x exp 



- / d*x E — 2-F + g( 7/ f D„ + m q )q + 8- 



4g2 - > ^ ^ ^ 1 1 32^ 2 



FF 



(C.l) 



For our purpose, let us assume only one heavy quark with mass toq. Then, 
under a chiral transformation 8^6 + 7r/2, the action will be invariant except for 
the following replacement: 



m Q -> m Q e l757r/2 = m Q ij 5 . 



(C.2) 



Appendix D. Heavy quark expansion 

Here, we rederive the heavy quark expansion for the lowest quark condensate^!. In 
the fixed point gauge, the heavy quark propagator can be written as an expansion 
in gluon operators. The term in the momentum space propagator that gives a gauge 
invariant term in the gauge invariant quark condensate operators is the following^!. 



iS G 2(p) 



where, 



[a, f3, p, a] = 



1 



[a, /3, p, a] + [a, p, /3, a] + [p, a, /3, a] 



1 



1 



1 



1 



$ + m ' a p 1 + m ^ + m ' p p 1 + m 1 a p 1 + m 



"7cr- 



Substituting, 



(qq) = - lim .Tr(T[q(x)q(0) > 



(D.l) 



(D.2) 



d A p 

WT 4 



(Tr[iS(p)}). 



(D.3) 



Substituting Eq. (|D. 1|) into the above formula, one finds the following leading term 
in the heavy quark mass expansion for the heavy quark condensate. 



(m h hh) 
(m h hij 5 h) 



.1(^G 2 ) + (A) 
12 it m h 

l{-GG) + 0(-L) 
8 7T m 2 h 



(D.4) 
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